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ABSTRACT 

Green Peas are a new class of young, emission line galaxies that were discovered by citizen volunteers 
in the Galaxy Zoo project. Their low stellar mass, low metallicity and very high star formation 
rates make Green Peas the nearby (z ~ 0.2) analogs of the Lyman-break Galaxies (LBGs) which 
account for the bulk of the star formation in the early universe {z ^ 2 — 5). They thus provide 
accessible laboratories in the nearby Universe for understanding star formation, supernova feedback, 
particle acceleration and magnetic field amplification in early galaxies. We report the first direct 
radio detection of Green Peas with low frequency GMRT observations and our stacking detection 
with archival VLA FIRST data. We show that the radio emission implies that these extremely young 
galaxies already have magnetic fields (> SO/xG) even larger than that of the Milky Way. This is at 
odds with the present understanding of magnetic field growth based on amplification of seed fields 
by dynamo action over a galaxy's lifetime. Our observations strongly favor models with pregalactic 
magnetic fields at levels. 

Subject headings: galaxies: magnetic fields — radio continuum — galaxies: starburst — cosmic rays 



1. INTRODUCTION 

A new cla ss of young em ission line galaxies, named 
Green Peas ([Cardamone et al. 2009i) , was recently dis- 
covered by volunt eers in the Galaxy Zoo project 
(jLintott et al.l [20081 ). They are particularly interesting 
because of their small size and their extraordinarily large 
[Oiii] equivalent width (up to ^1000 A; see Fig. [T|), which 
indicates a large population of hot young stars. They 
turn out to be low mass galaxies (M^ 10®-^ — 10^° M©) 
with high star formation rates (~ 10 M0yr~^). Hence 
they have some of the the highest specific star forma- 
tion ratefl seen in the local Universe. This indicates 
very short stellar mass doubling times (^double ^ lO^yrs, 
which is ^ 10~^ times the Hubble time), whereas typical 
star forming galaxies have stellar mass doubling times 
comparable to the Hubble time. Due to their small 
reddening (E(B — V) < 0.25) and luminous UV emission 
this class of star forming galaxies is similar to UV lu- 
minous galaxies observed at high redshift (z ~ 2 — 5). 
Their relatively low metaUicities (log[0/H] + 12 --8.7) 
and locations in under-dense environments suggest that 
Peas started building their stellar content at z « 0.2 by 
processes similar to those experienced by high redshift 
Lyman-a emitters and Lyman Break Galaxies (LBGs; 
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^ M /M of up to ~ 10~** yr-i 
ll2009l). Typical Peas h ave r 
ICardamone et al.l l[2009t) . 



were obtained bv lCardamone et al.l 
~ 20, n - r ~ 0.9 and z ~ 0.2 



see IGiavaliscd [2(M fo r a review). It has been sug- 
gested bv lAmorfn et al.l ([2010f ) that Peas may have sys- 
tematically larger [N/0] ratios, and th erefore potentially 
even l ower metallicity than reported bv lCardamone et al.l 
(l2009l) . 

Though unresolved by the SDSS, requiring typical an- 
gular sizes <1" which gives physical sizes <5 kpc, some 
Peas had serendipitous Hubble Space Telescope (HST) 
observations. Their patchy irregular appearance in re- 
solved HST s napshots possibly due to compact star form- 
ing regions (ICardamone et al.l [20091) . mimics the mor- 
phologies of high redshift galaxies. Given high star for- 
mation rates, Peas are expected to host a large num- 
ber of supernovae. Supernovae accelerate electrons to 
high energies, which may then emit synchrotron radia- 
tion in radio bands. We have performed a stacking ex- 
periment with arch ival Very Large Array (VLA) FIRST 
([Becker et al.lll995[ ) data at 1.4 GHz to demonstrate that 
their fiuxes are systematically different from those of 
usual starburst galaxies but are statistically akin to those 
of LBGs, which have systematically lower fi uxes than lo- 
cal st arbursts at a given star formation rate ([Carilli et al.l 
I2008D . We have followed up three Peas with deep ob- 
servations at 0.6 GHz using the Giant Metrewave Radio 
Telescope (GMRT). In this letter we establish and discuss 
the properties of Peas as a new class of sub-milliJansky 
sources in the GHz radio sky. We use theoretical consid- 
erations to show that these radio detections imply large 
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Fig. 1.— SDSS spectra of three Peas, targeted with GMRT. All 
spectra (in black) are de-redshifted and fitted (in red) with the 
Gas And Absorption Line Fitting (GANDALF; ISarzi ct al.][200g) 
that simultaneously fits both stellar spectra and emission lines to 
det ermine the metallicity, s tellar mass and star formation rates. 
See lCardamone et al.l II2009I ') for details of the procedure. 

magnetic fields for Peas. Tliis raises puzzling questions 
about the origin and evolution of magnetic fields in young 
galaxies. 

2. RADIO EMISSION IN GALAXIES 

Radio emission from most normafl galaxies is due to 
synchrotron radiation of relativistic electrons and free- 
free e mission of ionized hydrogen (Hii) regions (jCondonI 
|1992|) . Both of are produced by massive stars (M > 
8Mq) which end their bright but short lives (< 3x 10'' yr) 
in supernovae. During their lifetimes these stars ionize 
the surrounding Interstellar Medium (ISM) and create 
Hii regions. Shocks produced by supernovae accelerate a 
majority of the relativistic electrons in normal galaxies. 
These electrons have short lifetimes (< 10^ yr); hence 
radio emission in a normal galaxy traces recent star for- 
mation activity. 

The non thermal radio brightness of a typ ical star 
forming galaxy is given by lYun fc Carillil ()2002f ) as 



= 25 /„thJ^-° 
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where is luminosity distance in Mpc. Variations in 
the normalization of the nonthermal synchrotron emis- 
sion are traced by /nt h, which is ^ 1 for t he Milky Way 
and local starbursts (|Yun fc Carillil I2002D . With high 
star formation rates, young star forming galaxies like 



^ Here normal IICondonlll992l 'l refers to nearby galaxies in which 
the radio emission is n ot dominated by a cen tral super-massive 
blackholc. Scyfcrt Peas llUardamone et al]|2009l ) are excluded from 
this study and we only concentrate on Star Forming Peas. 
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Fig. 2. — Expected Signal (in mJy), 6a detection threshold (in 
mJy) and Signal to Noise Ratio (SNR) for a simulated stacking 
experiment with N brightest Peas. Note that the SNR rises and 
then falls with N, so stacking all available Peas is not optimal. We 
stack 2'' Peas in this work. 

LBGs and Peas are expected to have some nonthermal 
radio emission. However, to date no direct radio emis- 
sion has be en reported. In LBGs, stacking techniques 
(jCarilli et a l. 2008) have been used to overcome large dis- 
tances. With Peas, we have a local sample to search for 
their radio emission both by stacking existing VLA data 
(Section [3]) and by additional observations with GMRT 
(Sectiongl. 

3. VLA FIRST STACKING DETECTION 

Highest estimated fluxes for Peas fall jus t short of the 
VLA FIRST survey threshold of 1 mJy (|Becker et al.l 
|1995[ ) at 1.4 GHz, so they ne ed to be detected by stack- 
ing. See I White et al.l ()2007[ ) for extensive discussion of 
techniques for stacking FIRST images to detect sources 
below the survey threshold. Here we use a summation 
based approach, rather than a median based approach, 
as this gives higher Signal to Noise Ratio (SNR). This 
is because, the convergence of the sample mean to the 
population mean is unbiased and faster than the conver- 
gence of the sample medi an to the population median 
(jKennev fc Keepinel 119571) . Conventional wisdom sug- 
gests that summing up N images of faint objects produces 
a increase in SNR by a/TV This implies that stacking all 
available objects is the best strategy. However, this strat- 
egy can be bettered if some flux estimate is available. 
We estimate the expected fluxes of all 80 spcctro scopi- 
cally i dentifled star forming Peas from ICar damonc et al.l 
(|200a using Eq|T]with /nth = 1 in a flat ACDM cosmol- 
ogy with h = 0.71 and JIm = 0.3, with star formation 
rates derived from their optical spectra (See Figure |T]). 
In this situation it is best to stack the N galaxies with 
the highest expected fluxes. In Figure [5] we show that 
the SNR rises rapidly in the range TV ~ 1 — 15 and falls 
off slowly beyond > 30. We therefore stack A^ = 2^ 
(the nearest power of 2) of the brightest expected Peas. 

The images of these A^ regions were cut out and 
summed in pairs, with this process being repeatedly ap- 
plied to the resulting A^/2 maps until the flnal map (Fig. 

was obtained. This was done to avoid a sequential 
summation, which results in the addition of large floating 
point numbers with small ones in the latter steps, lead- 
ing to accumulation of rounding errors. In our method, 
the numbers being added are always comparable. The 
final image (Fig. [S]) has noise which is well described by 
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PLot file version 1 created 13-APR-2011 03:52:23 
BOTH: J082247+ GP SMOOTH. C0NVL.1 
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Center at RA 08 22 47.743 DEC 22 41 42.78 
Grey scale fiux range= 0.000 4.000 MiiiiJY/BEAM 
Cent peait fiux = 8.4766E-03 JY/BEAM 
Levs = 8.545E-04 * {2, 3, 4) 



Fig. 3.— VLA FIRST detection of stacked flux from 32 Peas. 
The displayed map is smoothed to a restoring beam size of 7". 
The unresolved 4(t source in the centre is due to the Peas. The 
only other significant detection, the resolved source to the right of 
the image, is due to the unrelated FIRST J1SS926. 7+151655. 

a Gaussian with a = 0.861 mJy. The stacked total flux 
is found to be 3.975 mJy instead of the expected 7.432 
mJy for /nth = 1, summed over all N Peas using Eq[T]for 
a standard starburst template. The probability of a sta- 
tistical fluctuation of this size and larger is ^ 9.2 x 10"'*. 
This indicates that the average radio fluxes of Peas are 
systematically suppressed, compared to usual starburst 
galaxies, by a fact or of /nth — 0.53. Thus, as in LBGs 
(jCarilli et al.ll2008h . Peas have a comparable but system- 
atically lower flux when compared to local starbursts. 

4. GMRT OBSERVATIONS AND RESULTS 

The result from stacking demands confirmation at 
other frequencies, especially in terms of direct radio de- 
tection of Peas. Hence, we targeted the 3 most promising 
candidates, with expected fluxes at the mJy level, using 
deep GMRT observations. We detect two of them as 
unresolved point sources and put an upper limit on the 
third. These results show that not all Peas may be char- 
acterized by a single /nth- Observations with ~ 3 hours 
of on-source time for each Pea, were carried out at 617 
MHz (bandwidth 32MHz, split into 512 channels). Visi- 
bility data have been analyz ed using Astro nomical Image 
Processing System (MPS) (IGreisenlll99nl) . 

Initial calibrations of the data were done using various 
flux calibrators and phase calibrators (See Tabled]). The 
flux calibrators were observed at the beginning and the 
end of the respective observations in order to fix the an- 
tennae gains. The phase calibrators, chosen within 15° 
of the source position on the sky, were observed every 
half an hour to calibrate the phase drift due to iono- 
spheric effects. The data were edited for the removal of 
bad antennae, baselines and any radio frequency inter- 
ference (RFI). The data were binned into seven channels 
to maximize SNR without compromising on chromatic 



aberration. Target visibilities were extracted from gain 
and bandpass calibrated data. The task IMAGR was 
used on this data to obtain a preliminary dirty map. 
Cleaning was done to remove the beam pattern for all 
3(7 and brighter sources. Repeated self calibration and 
cleaning was done till it yielded consistent solutions. The 
source flux was then extracted with the task JMFIT from 
the final map, with a typical resolution of ~ 6" . Table 
[1] lists sources, calibrators, fluxes and implied magnetic 
fields. 

5. ELECTRON DIFFUSION 

Relativistic electrons produce bulk of the GHz radio 
emission in star forming galaxies. They must remain 
confined long enough in the galaxy's magnetic field to 
significantly contribute to the radio emission. We use 
the criteria derived here, to constrain the characteris- 
tic magnetic fields in Section [T) These electrons trav- 
eling at speeds close to that of light get scattered by 
the galaxy's magnetic field. Their diffusion coefficient is 
determined by t he structure of the m agnetic field. We 
assume foUowing lGinzburg &: Svrovatsk ii (196^) that the 
field structure is comprise d of domains of rad ius Iq (typi- 
cally 0.3 kpc according to lCalvez et al.l ()2010f )) with ran- 
domly oriented magnetic field B in each of these domains. 
For the range of feasible electron energies, the Larmor 
radius is smaller than the magnetic coherence length, 
thereby providing an energy independent diffusion co- 
efficient D = Iqc/3. Therefore the diffusion timescale 
for relativi stic electrons to escape from a typical Pea of 
size 3 kpc (jCardamone et al.]|2009() . is 

This can be compared to the timescales of radiative 
processes to gauge their relative importance. 

6. INVERSE COMPTON SCATTERING 

Electrons may lose their energy before diffusing out of 
a galaxy via Inverse Compton (IC) scattering. An ex- 
planation of t he systemati c ally lo wer fiuxes of the LBGs 
suggested bv ICarilli et al.l (|2008D . is IC cooling of rela- 
tivistic electrons scattering off Cosmic Microwave Back- 
ground (CMB) photons. If this is indeed the cause, 
the effect should be negligible for Peas as the CMB ra- 
diation density in the Universe has fallen off consider- 
ably in the epoch between LBGs (z ~ 2 — 5) and Peas 
(z ~ 0.2). Electrons can lose their energy through IC. IC 
has the same dependence on the photon energy density 
as synchrotron has on the m agnetic field energy density 
()Rvbicki fc Lightman|[l979l ). This implies the IC energy 
loss timescale; 



~1 29 X 10^ vr [ — ^ ( ^] '^'"^ 

^M^lO-i^ergcm-^y VCHz^ 

where C/rad is the energy density in radiation or photons. 
For IC losses to be significant, t\c has to be comparable 
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TABLE 1 
GMRT Observations of Peas 



Pea 

(SDSS Object) 


z 


Flux Cal 
(3C Source) 


Phase Cal 
(J2000) 


(mJy) 


AtG 


Bcq 


J082247.66+224144.1 
J074936. 76+333716.3 
J142405. 73+421646.3 


0.216 
0.273 
0.185 


3C147, 3C286 

3C147 

3C286 


J0741+312 
J0842+185 
J1416+347 


1.20 + 0.31 
1.11 + 0.11 
< 0.45 


81 + 14 

82 + 5 
< 49 


50 + 4 
54 + 1 
< 35 



Note. — Details of three Peas observed at 617 MHz: Bj^ and Bcq are the magnetic fields 
derived using diffusion and equipartition arguments. Uncertainties are Icr statistical standard 
errors. The entries for the last object represent the 3cr upper confidence limits. 



to to- This condition gives the required i/rad as 
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The energy density of the CMB photons is given by 



CMB 



= aT'' ~ 4 X 10"^^(1 + zf erg cm 
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This gives the time scale for IC losses against CMB pho- 
tons as 



^CMB 



E 



\dE/dt\ic 



(6) 



.2.55xl0«yr {I + z)-' (^)"'^'. 



Comparing this to we get the redshift, at which the 
electrons lose a significant fraction of their energy to 
CMB photons, as 



1 + z ~ 9.54 X 



R 
3kpc 



-2/3 



-1/6 



0.3kpc 



1/3 



(7) 



This may be a significant effect for very high redshift 
(z ~ 8.5) galaxies. Thus CMB IC cooling of rclativistic 
electrons may have a role in explaining the systematic 
lowering of LBG radio fluxes. However it cannot have 
any significant role in the case of Peas which are at z ~ 

0.2^ 

iCondon et al.l (|1991| ) suggested that radio emission 
from compact starbursts in ultra-luminous infrared 
galaxies may be suppressed due to IC losses against the 
radiation energy density produced by its own young stel- 
lar population. Comparing the IC loss timcscalc to the 
electron diffusion timcscalc, we have obtained the re- 
quired energy density in Eq. |4l The typical value of ?7rad, 
required for IC losses to significantly suppress the syn- 
chrotron flux at G Hz frequencies, is much larger than the 
Milky Way value (jCondonl 11993) of IQ-^^eig cmT^ but 
less than t he higher value s of up to 10~*erg cm~'^ which 



are seen bv ICondon et al.l (|1991| ) in compact starbursts. 
So it is possible that some Peas may have suppressed 
radio fluxes as a result of IC losses against their own ra- 
diation density leading to a wide range of observed values 

for /nth- 

One out of three Peas observed by GMRT showed no 
detectable flux (Table [1]). We therefore estimate [/lad in 



this galaxy from its SDSS spectra, applying appropriate 
bolomctric correction for a starburst galaxy. We find an 
t^rad ~ 4 X 10~^^erg cm"'^ which rules out IC cooling of 
electrons as the cause for suppression of radio luminosity. 
Hence the upper limit of flux from this object allows 
us to limit the magnetic field using the diffusion and 
equipartition arguments. 

7. MAGNETIC FIELD: DIFFUSION ARGUMENT 

To explain the observed radio emission from Peas as 
synchrotron loss from electrons, the magnetic field must 
be large so that electrons lose enough energy before they 
diffuse out. This can be used to constrain magnetic 
fields in Peas. Since CMB IC cannot account for the 
systematic lowering of Pea radi o fluxes we explore an- 
other explanation suggested by iCarilli et al.l (|2008[ ) for 
the suppression of LBG fluxes. Cosmic Rays may dif- 
fuse easily from systematically smaller galaxies in the 
early Universe before they can lose their energy via syn- 
chrotr on emission. Since Peas are of sizes comparable to 
LBGs (jCardamone et al.] 120091 ) this effect should be ob- 
servable in Peas. Here we use the observed radio emis- 
sion to derive the ma gnetic field of a Pea, by compa ring 
the electron diffusion (jGinzburg fc Svr ovatskii"1964') and 
synch rotron energy loss timescales (jRybicki & Lightma^ 
Il979( ). Thus we get a magnetic field required to explain 
the observed radio flux; 



Bd ^39 nG 



( R 



\3kpc 



Mgh^) 



-4/3 

/nth 
ll5" 



0.3kpc 

2/3 



2/3 



(8) 



where the magn etic coherence leng th Iq has a typical 
value of 0.3 kpc (jCalvez et aLllMO) . Using /nth ^ 0.53 
at Usy-a = 1.4 G Hz (from stacking exper iment) and char- 
acteristic size (jCardamone et aLl l2009f ) i? of 3 kpc, we 
obtain a typical magnetic field oi B k, 36^G. Thus the 
electron diffusion argument provides us with a direct es- 
timate of the characteristic magnetic field in Peas. 

8. MAGNETIC FIELD: EQUIPARTITION ARGUMENT 

We show here that the previous estimate of magnetic 
field does not lead to absurd energy requirements. An 
independent estimate of magnetic field can be obtained 
from energy minimization considerations. The so called 
equipartition magnetic field can be estimated by mini- 
mizing total energy invoked in rclativistic particles and 
magnetic fields, in order to explain the observed flux. 
The magneti c field corr e spond ing to the minimum energy 
condition of iBurbidg^ (|1956f ) in a synchrotron plasma 
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can be deriv ed from radio obse r vation s. We follow the 
derivation of iFitt fc Alexander! (|1993D with the typical 
value of the cosmic ray prot on/electron energy ratio 

jical 
for 

the radio spectrum, to derive the equipartition magnetic 
field as 



value of the cosmic ray prot on/electron energy rat 
(jGinzburg fc Syrovatsknlll964D k = 100. and a typ ic 
spectral index of a 3/4 (|Fitt fc Alexander! Il993f ) f 



/ c \ 2/7 



-4/7 



VGHz/ 



(9) 



If the actual magnetic field is different from this value, 
the energy required to explain the synchro tron emission 
goes up very rapidly. It has been argued bv lDurid (|199Clf ) 
that the real magnetic field is limited to differ at most by 
an order of magnitude above and below the equipartition 
value derived in this manner. From the stacking experi- 
ment the mean flux of the Peas is ^ 0.12 mJy. While 
they are unresolved at VLA FIRST re solution, serendip- 
itous HST observations (jCardamone e t al. 2003) provide 
their characteristic size as 3 kpc (where 1" is ^ 5 kpc at 
their typical redshift). 

This fixes the characteristic equipartition magnetic 
field of Peas at 39/iG. The surprising agreement be- 
tween the estimates of Pea magnetic fields derived from 
two independent methods, points to a characteristic mag- 
netic field of > 30^G which is comparable to or even 
larger than th e average Milky Way value of S w 5/iG 
(|Gondonlll99"2h . 

9. DISCUSSION 

We report the first radio detection of Peas and show 
that this discovery implies large magnetic fields in Peas, 
using reasonable assumptions about Cosmic Ray diffu- 
sion and total energy considerations. Given that the 
Peas are very young galaxies (around one hundredth the 
« 10^" yr age of the Milky Way) the discovery of radio 
emission and the implied fiG magnetic fields is a striking 
result. Present day magnetic fields are thought to be the 
result of amplificat ion of seed fiel ds (~ 10^^° — 10~^^G) 
by dynamo action (IWidrowll2002D over the galaxy's life- 
time. These models for the growt h of magnetic fields 
have e-folding tim es jW idrow 2002") of sa 1 0^ - 1 0^ yr 
They produce 20 (jKukrud fc Zweibe1[200l to 50 (|Ree^ 



l2006f ) e-folds of an efficient dynamo during the lifetime 
of the Milky Way. Hence, they can only play a small role 
in amplifying the pregalactic (jKulsrud fc Zweib"eill2008D 
magnetic fields within the age of a Pea. It has been sug- 
gested bv lDalv fc Loebl (|1990( ) that galaxy- wide fiG level 
magnetic fields may be generated by magnetized plasma 
produced in jets from a central compact object. Even 
this is unlikely as we are considering only the Star Form- 
ing Peas and not the Se yfert Peas. Our resu l ts str ongly 
favor the suggestion from lKulsrud fc Zweibell (|2008| ) that 
seed fields were amplified significantly by turbulence (up 
to level pregalactic fields) as protogalaxies and simi- 
lar substructures formed. 

As Peas are close analogs to LBGs it seems likely 
that young primeval galaxies could have had strong 
enou gh magnetic fields to se ed the intcrgalactic magnetic 
field iKronberg et al.l [l999h . It has been suggested by 
iZeldovich et al.l (|1983r ) that magnetic fields may play a 
crucial role in star formation. Observed magnetic fields 
in ne arby galaxies scale slowly with the SFR (jValled 
|1994[ ) but it is likely that magn etic fields str ongly af- 
fect the star formation efficiency (|Totanilll999( ). In such 
a scenario, presence of strong pregalactic magnetic fields 
may be the reason why Peas have started to produce 
stars at such a large rate. Most theoretical and compu- 
tational studies of reionization er a assume lack of dynam - 
ically significant mag netic fields (|Loeb fc BarkanaH200lD . 
during formation of first stars. Our results may require 
a rethinking of this assumption, in presence of fiG level 
pregalactic magnetic fields. A detailed study of processes 
at work in Peas in the nearby Universe, will help us un- 
derstand magnetic field amplification in early galaxies, 
its role in star formation, supernova feedback and cos- 
mic ray acceleration. 
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